The rate of stellar tidal disruption flares from SDSS data 
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Abstract We have searched for flares due to the tidal disruption of stars by supermassive black holes in archival 
Sloan Digital Sky Survey (SDSS) multi-epoch imaging data. Our pipeline takes advantage of the excellent as- 
trometry of SDSS to separate nuclear flares from supernovae. The 10 year baseline and the high cadence of the 
observations facilitate a clear-cut identification of variable active galactic nuclei. We found 186 nuclear flares, 
of which two are strong stellar tidal disruption flare (TDF) candidates. To compute the rate of these events, 
we simulated our entire pipeline to obtain the efficiency of detection for a given light curve. We compute a 
model-independent upper limit to the TDF rate of N < 3 x 1(T 4 yr~'galaxy~' (90% CL). Using a simple model 
to extrapolate the observed light curve forward and backward in time, we find our best-estimate of the rate: 
N = 3 + J 3 x lO-'yr^galaxy" 1 . 



1 Introduction 

A star that comes too close to a massive black hole is torn 
apart by tidal gravity forces, yielding a stellar tidal dis- 
ruption flare (TDF). How many stars per galaxy per year 
suffer this fate is currently not well constrained (see [ 1 1 for 
a recent review). Based on the nuclear density profiles of 
nearby elliptical galaxies, Wang and Merritt [2] compute 
a rate of ~ 10~ 4 yr , with 1 dex scatter between different 
galaxies of similar mass. 

Recently, two TDFs were found in the Sloan Digital 
Sky Survey (SDSS) multi-epoch imaging data [3, here- 
after paper I]. The systematic nature of this search allows 
for a relatively straightforward computation of the rate of 
these events. In this proceeding, we shall first summarize 
the search of paper I, followed by a derivation of the TDF 
rate. Full details on rate analysis will be given in a forth- 
coming publication. 



1.1 Summary of SDSS nuclear flare search 

The search for TDFs (paper I) was conducted in Stripe 82. 
Nuclear flares in galaxies are found using two steps: (i) a 
series of catalog cuts to select flaring galaxies, followed 
by (ii) careful difference imaging to measure the angular 
distance between the flare and the host. The catalog cuts 
selected galaxies with a flux increase of 10% or more, de- 
tected at the 7<r-level. A sample of 186 nuclear flares was 
selected based on the distance between the center of the 
host and the center of the flare in the difference image 
(d < 0.2"). After removing galaxies that fall inside the 
photometric QSO locus and removing galaxies with addi- 
tional variability, two flares remained: TDE1 and TDE2. 
Additional analysis and follow-up observations show that 
these flares are best explained as stellar tidal disruption 
events (paper I). 
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2 Analysis 

The number of detected flares in a variability survey that 
targets galaxies is given by 

A%jf = t ^ e,^,- (1) 

i 

where N, and e, are the flare rate and detection efficiency 
for the /th monitored galaxy, and r is the survey time. The 
rate of TDFs is expected to depend only weakly on black 
hole mass [4|. Furthermore, as we shall show below, our 
search is sensitive to a relatively narrow range of black hole 
masses. We can therefore simplify Eq.[T] using N t — N, a 
galaxy-independent rate, to find 



Here we defined we mean overall efficiency e = N 1 2f 
For the TDF search in Stripe 82 we set t = 7.6 yr, starting 
in the year 2000, Ntdf = 2, and N g - dl = 1.6 x 10 6 , i.e., 
all galaxies with a photometric redshift that are outside the 
QSO locus. Finding N thus boils down to computing the 
mean overall efficiency, e. 

As discussed in sec. |1.1| the detection pipeline of pa- 
per I consists of two stages: the catalog cuts and the dif- 
ference imaging. Since the catalog cuts are applied to the 
Petrosian flux [6| of the galaxy, computing the probabil- 
ity that a simulate light curve passes these cuts is trivial: 
one simply adds the flare flux to the galaxy flux and re- 
runs the catalog cuts. To estimate the detection probability 
of the difference imaging pipeline we selected 1400 ran- 
dom galaxies in uniform magnitude bins and inserted point 
sources at the center of their images. The detection prob- 
ability as a function of flare and host magnitude follows 
from the number of detected point sources in each magni- 
tude bin. 

With the detection probabilities measured, we now pro- 
ceed to compute the overall efficiency (e in Eq. pp. For 
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Figure 1. Sensitivity of our search (i.e., the denominator of Eq. 
[2] N^ite), in bins of absolute magnitude of the host (thin and 
thick lines for TDE1 and TDE2, respectively). The mean black 
hole mass in each bin is indicated on the upper axis. We also 
show the parent galaxy sample by setting e = 1 (thin black line). 
Our analysis is most sensitive to galaxies with M r = -21.5 or 
M BH ~ 5 x 1O 7 M . 

a given galaxy we first draw a location for the start of 
the flare for a uniform distribution. We then add the flare 
flux to the Petrosian flux and check if this galaxy would 
pass our catalog cuts. In the final step, we simulate the 
detection of this flare (for at least three nights in the u, 
g, and r bands) using the probability of detection for the 
flare and host magnitude. After repeating this process for a 
sufficiently large sample of galaxies, the overall efficiency 
follows from the number of flares detected by the model 
pipeline over the number of simulated flares. Because the 
flares are inserted into the observed galaxy light curves, 
this method fully takes into account the inhomogeneous 
cadence and varying data quality of Stripe 82. 

The detection efficiency will obviously depend on the 
flare's luminosity and duration. Because the SDSS data of 
the two TDFs of paper I does not completely cover the time 
they are detectable (i.e., above the flux limit), we have to 
extrapolate the observed light curve forward and backward 
in time. Since our optical observations probe the Rayleigh- 
Jeans regime of the SED, we adopt L TDF oc f~ 5 / 12 [5|. We 
fit the observed light curve to F oc (t - toT 5 ^ 12 to find the 
time of disruption, tp. To fix the normalization of the light 
curve, we shall assume that the luminosity of the flare is 
proportional to the Eddington luminosity, Ltdf k Mbh- An 
estimate of the black hole mass (Mbh) is obtained from the 
galaxy luminosity using the black hole-bulge mass relation 
Q and the dynamical mass-to-light ratio from the funda- 
mental plane, yielding Mbh k L^i e' w ^ -^buige the bulge 
luminosity. In summary, the expression for the model light 
curve of a flare in the z'-th galaxy is 

Mi(t) = MtdE1,2(0 +1-3 (M ; , bu i ge - M TD El,2,bulge) (3) 

with Mtdei,2(0 the (extrapolated) light curves of TDE1,2. 
We use the galaxy photometric redshift 1 8 1 to convert be- 
tween apparent and absolute magnitudes. 

The light curves we simulate are 300 days long, ex- 
trapolated back in time to five days after the disruption, 
to — (5 + to) (we do not extrapolate back further because at 
very early time the power-law scaling is not appropriate). 



Table 1. Measured rate of optical stellar tidal disruption flares 

Light curve Rate 

(Ftdf « t p ) (galaxy -1 yr~') 

Upper limit Observed < 3 x 10" 4 

Disk model P = ~j2 3 -3 x 10 " 5 

The final efficiency is quite insensitive to the light curve 
length or to because we required at least three detections 
of the flare within the same season. 

Finally, we also compute a model-independent upper 
limit to the TDF rate by using only the observed fraction 
of the light curve of TDE1,2. 



3 Results 

Using the model light curves presented in the previous sec- 
tion (Eq. [3j, we injected flares scaled to either TDE1 or 
TDE2 to find the overall efficiency for detecting these flares. 
We use the average of these two overall efficiencies to ob- 
tain the TDF rate (Eq.|2]i; the results are shown in Table[T] 
The uncertainties on the rate are given by the 90% confi- 
dence level (CL) for Poisson statistics; for the upper limit 
on the rate we use Njdf < 5.3, the 90% CL upper limit if 
two events are detected. 

In Fig. [T]we show the sensitivity of our pipeline as a 
function of host luminosity. The sensitivity to TDFs from 
faint M r < -20 galaxies decreases faster than their ob- 
served number density, which is due to the scaling of the 
flare luminosity with black hole mass. Our search is most 
sensitive to galaxies in the range -21 < M ( - < -22.5. For 
each host we can estimate the central black hole mass from 
the bulge luminosity using the black hole mass-bulge mass 
relation, calibrated for the SDSS r-band 0. We find that 
our analysis is most sensitive to black holes in the mass 
range 2 x 10 7 M to 3 x 10 8 M . 

The TDF rate we derive with the simplified model light 
curve is below some of the theoretical predictions that have 
been published J2), but consistent with the rate of soft X- 
ray flares from quiescent galaxies [ 101 . Our strictly obser- 
vational upper bound is not strong enough to rule out any 
current predictions. 



References 

1. T. Alexander, arXiv e-prints (2012) 1210.8582 

2. J. Wang, D. Merritt, ApJ 600, 149 (2004) 

3. S. van Velzen, G.R. Farrar, S. Gezari, N. Morrell, 

D. Zaritsky, L. Ostman, M. Smith, J. Gelfand, A.J. 
Drake, ApJ 741, 73 (2011) 

4. J. Magonian, S. Tremaine, MNRAS 309, 447 (1999) 

5. G. Lodato, E.M. Rossi, MNRAS 410, 359 (2011) 

6. C. Stoughton et al., AJ 123, 485 (2002) 

7. N. Haring, H. Rix, ApJ 604, L89 (2004) 

8. H. Oyaizu, M. Lima, C.E. Cunha, H. Lin, J. Frieman, 

E. S. Sheldon, ApJ 674, 768 (2008) 

9. E. Tundo, M. Bernardi, J.B. Hyde, R.K. Sheth, 
A. Pizzella, ApJ 663, 53 (2007) 

10. J.L. Donley, W.N. Brandt, M. Eracleous, T. Boiler, 
AJ 124, 1308 (2002) 



